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ABSTRACT 

An investigation has been made of the feasibility of making wind 
velocity measurements from spare by monitoring the apparent change in 
the refractive index of the atmosphere induced by motion of the air. The 
physical principle is the same as that resulting in the phase changes mea- 
sured in the Fizeau experiment. It is proposed that this phase change 
could be measured using, for example, a three cornered arrangement of 
satellite borne source and reflectors, around which two laser beams pro- 
pagate in opposite directions. It is shown that even though the velocity 
of the satellites is much larger than the wind velocity, factors such as 
change in satellite position and Doppler shifts can be taken into account 
in a reasonable manner and the Fizeau phase measured. This phase mea- 
surement yields an average wind velocity along the ray path through the 
atmosphere. The method requires neither high accuracy for satellite 
position or velocity, nor precise knowledge of the rafractive index or its 
gradient in the atmosphere. However, the method intrinsically yields 
wind velocity integrated along the ray path; hence to obtain higher spa- 
tial resolution, inversion techniques will be required. This paper addresses 
the general principle of the technique and presents a particular system 
configuration as an example, to show that wind measurements are 
possible. 
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A PROPOSED METHOD FOR WIND VELOCITY MEASUREMENT 


FROM SPACE 

INTRODUCTION 

The present report is a study of the potential for measuring wind velocity using the pnase 
change induced in two coherent laser beams traversing the atmosphere in opposite directions. It 
will be shown that a coherent measuremen "■ feasible from space and can be accomplished in spite 
of the motion of the satellites. At a first glance this seems to be a formidable task. Satellites are 
moving at velocities of the order of 10km/ sec, while wind velocities arr of the order of 1 
to 10m/ sec, a difference of several orders of magnitude. In addition, the position of the satellites is 
not available to within an accuracy of an optical wavelength, which would seem to be important for 
a coherent experiment. However, it is shown below that these problems can be adequately resolved. 
This method does not require extraordinary accuracy of satellite position and velocity, in fact, 
resolution of the order of l-10m is sufficient. Also, one does not need a detailed knowledge of the 
atmospheric refractivity and its gradient, beyond that available from existing models and data. On 
the other hand, we have assumed the availability in space of laser systems for remote sensing pur- 
poses. While this is still a thing of the future, there is a reasonable prospect that such systems will 
exist towards the end of the century. 

In order to put the present study in the right perspective, it is worthwhile to quote from Gor- 
don Little* : ", . . In general, the development of a remote sensing concept can be seen to follow a 
logical sequence. In Step A, the concept is identified, and preliminary first order estimates made of 
its feasibility. In Step B, the potential capabilities and limitations of the concept are analyzed theo- 
retically in considerable detail. If the concept still appears attractive, the development of research 
equipment for the experimental evaluation of these capabilities and limitations takes place in Step 
Assuming that this quantitative experimental evaluation of the concept is successful, the next 
stage (Step D) is to build a development model (as opposed to a research model) which is thought 
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of a* a prototype of an operational unit which is to be capable of being used in the field by the 
research workers or technicians other than the original research group. If the concept continues to 
show promise. Step E involves working with industry to obtain commercially built units for field 
evaluation. Once this stage has been successfully completed, the final stage (Step F) requires that 
fully evaluated commercial units be routinely available for procurement.” Obviously the present 
report is mainly concerned with Step A. At this stage we close our eyes to budgetary requirements, 
and important engineering problems such as detectability of signals by means of existing or pro- 
jected optical instrumentation, the problem of tracking, and probably many others. 

The method to be presented below possesses many features which make it attractive as com- 
plementary to other systems when they operate under adverse conditions. For example, some 
information about winds on be obtained by monitoring cloud motion. In contradistinction, the 
present method is suitable for clear skies. Other microwave radar and radiometry methods' corre- 
late wind measurements and sea state. The present method will work equally well over sea and land 
and at arbitrary elevations. Unlike atmospheric radar or lidar, the present method does not depend 
on backscattering Ii ' a irregularities and particles. Since it is based on forward propagation, a 
transparent line of sight will yiekl the best results. Other line of sight methods exist, see for exam- 
ple Ishtmaru^ , relying on the presence of atmospheric inhomogeneities. These methods are particu- 
larly useful for measuring transverse winds, while the present method yields the wind component 
along the beam. It is also notable that the present method does not require exceptional imaging 
qualities, which would tend to make the op. ical equipment costly and heavy. 

In the sections to follow, the theory of the Fizeau expet iment is briefly reviewed and its appli- 
cability to ti e atmosphere considered. Then, typical configurations applicable to measuring winds 
are discussed. The Doppler effect is discussed in detail. This is relevant to the present method be- 
cause ot the change of frequency anti direction of propagation, occunng during reflection from 
moving satellites. Next, an outline of the signal processing necessary to effect a measurement is 
presented, showing that in principle the Fizeau phase offsets are measurable in the presence of the 



various Doppler effects. Finally a brief discussion of the coherence problem is given. Clearly, if 
the present concept survives scrutiny by the scientific community, this and many other aspects of 
the system will have to be discussed in greater detail. 
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THE FIZEAU EXPERIMENT AND ITS APPLICABILITY TO ATMOSPHERIC 
MEASUREMENTS 

The present method for wind velocity measurement is based on the measurement of the 
Fresnel convection coefficient in the Fizeau experiment 4 . One possible variant of the Fizeau exper- 
iment is depicted in Fig. 1 . This is an interferometer in which two light beams, emanating from a 
common source, traverse a moving fluid in opposite directions. Experimentally one finds that the 
emerging waves differ in phase by an amount which is proportional to the velocity. The analysis of 
the Fizeau experiment is based on the Lorentz transformation for frequency and propagation vec- 
tor. which to the first order in v/c takes the form 5 
k ' = k + u> v/c 2 . 

u>' = u> + k • v , ( 1) 

<"«w 

where v is the velocity of the fluid as observed in the laboratory frame of reference, c = 3 • 10® 

■*»» 

nVsec is the speed of light in free space; k and oo arc the propagation vector and (angular) fre- 
quency, respectively, in the comoving frame of reference where the medium is observed to be at 
rest. The primes denote quantities measured in the laboratory reference frame, quantities in the 
comoving frame have no primes. In the comoving frame the retractive index is defined by 



and n is assumed to be independent of frequency ( i.e. . nondispersive). Hence in the laboratory 
frame one observes, to the first order in v/c: 



(3) 


and when k and v are codirectional, this becomes 

n’ - n ( n* I ) v c, (4) 

wlieie it* I is the Fresnel convection coefficient Fins shows that the retractive index is at lei ted 


by the velocity of the moving fluid It is this etfect winch provides a method lor direct measuiement 


4 


of wind velocity in the atmosphere. The phase change accumulated by an electromagnetic wave 
along a trajectory (ray) is given by 



where the integration is along the ray from Pj to P 2 . Referring to Fig. I , it is clear that the total 
phase difference between the rays is 

4trLv , 

top = — — (n J - 1). (6) 

Xc 

For tenuous media like the atmosphere it is convenient to define a refractivity N by N * (n 0 • 
10 6 , hence (6) can be approximated by 

= 8rrLvN lO^AXc). (7) 

In an inhomogeneous medium this would be replaced by 


A<j> - 


SfflQ - 6 

Xc 


*^path 


vNdL . 


(8) 


where v is the component of the velocity parallel to the ray path. An average value for vNL for a 
given atmosphere may be defined by equating (7) and (8). 

To get an idea of the numbers involved, consider light at X = 0.5ft and a medium having n = 1.3 
(e.g.. water). Taking v= 10 nysec and L = 1 m in (6) yields: Aip = 0.58 rad = 33°. Such a value 
is easily measurable in an interference experiment by observing the shift of the fringe pattern rela- 
tive to its position for v = 0. In the case of the atmosphere, let 10 m/ sec — 20 mph serve as a 
typical value for the wind velocity. The atmospheric refractivity is on the order of N = 300 (e.g., 
see Bean and Dutton 6 ). Consider a source in the IR band, with X = ]Qft. In this region strong 
stable CO-, lasers are available and a window exists for which the clear atmosphere is practically 
lossless. Setting the path length at L = 100 km yields in (7) the result - 2.5 rad = 144°. This 
is an easily measurable phase. If we keep distances on a scale of L - 100 km and the wind velocity 
reduces to about 2 mph, then wc still have A>p - 14.4° which is not difficult to measure. However. 
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if we now choose a source in the microwave region for example with X * IQ" 2 m then Ap will be 
smaller by a factor of 10 3 . This would seem to completely rule out the use of microwaves for the 


terrestrial atmosphere. However, in planetary atmospheres one may encounter distances, velocities 
and refraction indices for which the use of microwaves may be of interest. 
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SATELLITE CONFIGURATION 

To illustrate that a measurement of wind velocity might be possible using the Fizeau effect, we 
will consider a hypothetical satellite configuration. In principle two rays are needed, one traveling 
through the atmosphere downwind, and the other traversing the atmosphere in the upwind direc- 
tion. Fig. 2 shows a potentially suitable satellite configuration. In this example M is a master sta- 
tion on which the laser source is located and on which the processing of the returned signals takes 
place. The slave satellites Sj and S 2 serve only to reflect the rays in the desired direction. The ray 
paths are marked 1, 1'. 1" and 2, 2', 2" in an obvious way. Configurations are also possible in 
which the earth (e.g., the ocean) is one of the reflectors. But the problem of maintaining beam co- 
herence becomes critical for these configurations. Systems in which the master station or a reflec- 
tor is located on the earth’s surface, or on board an aircraft are also potentially feasible. The 
configuration illustrated in Fig. 2 will serve here to describe the principle. 

Assume for the moment that all elements of Fig. 2 are at rest except for the atmosphere. Then 
the basic configuration consists of two laser beams which traverse the paths of equal length 1.1'. 

1" and 2, 2', 2 " in opposite directions. In principle, this geometry is identical to the Fi/cau experi- 
ment. (Fig. I) each beam will have experienced a phase shift which is due to two terms: the round 
trip distance and the contribution due to the motion of the medium. If the round trip distances 
are identical, then the only phase difference will be due to the motion of the medium, and since the 
beams traverse the moving medium in opposite directions with respect to the wind direction, the 
contributions to phase difference due to the Fizeau effect (8) will add. The net phase difference 
between the beams will then be given by (8) and can be measured by appropriately mixing the two 
beams at the master station. It is important to note that a knowledge of the total path length 
around the triangle is not necessary. Only the ray path through the medium and the value of N on 
it are needed. The necessary accuracy for these parameters is independent of wavelength, i.e.. in 
(7) it is the relative errors AN/N. Al./L that allect the result. To gel .111 Hle.i of * lie path length I. 
through the atmosphere consider its effective height to be h = 5 km and let the radius of the earth 





He a * 6360 km, Then from Fig. 3 L * 2atgo and cot a * a/(a+h). Substituting for h and a 
yield* L as $00 km. Thus, if N and L are each known to within say 1%, e.g.. 5 km in a 500 km 
path and 3 N units for N * 300, then the error in computing the wind velocity will be on the order 
of 27< . Reasonable models for N and its gradient dN/dh are available*, and satellite position . 
available to an accuracy of meters or better. 

The preceding arguments apply only if all elements of the system are at rest. If the satellites 
M, Sj and are allowed to move relative to each other and to the earth, as will be true in practice, 
then several important problems arise. First of all, changes will occur in the frequency of the laser 
pulses when they are reflected from the moving platforms Sj and S,. It will be shown that there is 
a net Doppler frequency shift between beams which travel around the triangle < M S| ) in opposite 

directions. This net effect gives rise to a nonvanishing phase difference between the two returning 
signals. It will be shown below how this effect can be taken into account and a meaningful mea- 
surement of the li/eau phase change made. A second problem encountered when the system moves 
is that the path length need not be the same for beams which traverse the triangle in opposite direc- 
tions. The upwind and the downwind beams(Pig. 2) encounter the satellites S ( and S-. at diflcicul 
times lor pulses leaving M simultaneously. If the satellites are moving, this means they are encoun- 
tered at different positions. This change in position introduces a net path length difference between 
the two beams. This is a real problem, because a path difference on the c 'er ol one wavelength 
would completely obliterate the phase difference due to the hi/cau effect However, it will be 
show n below that all motional effects can be properly taken into account in the computation ol Un- 
wind velocity. I'hc fact that the laser beams become nonplanar along the path because of beam 
spread and that the satellites do not move on straight lines i due to orbital motion) w ill also be dis- 
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cussed below 



ANALYSIS Oh SA TELLITE MOTIONAL EFFECTS 


Generally speaking, relative motion of the components of a system introduces frequency, 
direction and amplitude shifts, in addition to affecting the positions of the various parts of the sys- 
tem. These manifests! ions of motion can influence the phase of the returning rays. Consequently 
a procedure must be devised for deriving the Fizeau phase change ( 8 ) in the presence of the spurious 
phase factors. In a proper relativistic treatment of electromagnetic plane waves in moving systems* 
there appears an amplitude effect of the first order in the velocity. This does not affect the phase 
and therefore, for purposes of this discussion, it can be neglecteJ. Doppler shifts in wavelength and 
frequency are first order effects in v/c and since (7). ( 8 ) are already of the first order in v/c. these 
effects are of second order importance in the computation of the Fizeau effect. However, they 
have to be taken into account. For example, since we will be adding the returning waves in order 
to measure the Fizeau effect, any net frequency difference will cause the waves to beat in and out 
of phase, complicating the measurement of A<> Consequently a careful relativistic treatment of the 
Doppler frequency shifts will be presented below to properly assess these effects 

Keiativistic Doppler effect: Let two inertial systems of reference move at a relative velocity v . 

MM* “■*» 

The ‘laboratory’ system of reference x', t' is now attached to satellite M. The ‘comoving’ system 
x . t refers now to any other part of the configuration. Gravitational effects are neglected. For the 
special case where the origins x' = 0 . x = 0 coincide at t ® t' - 0 , special relativity prescribes 
the transformation 7 •* 

x’ = II • (x + vt) , x = U • (x' vt') 

t' = 7 (t + v • x/e J ) . t = >(t' v • x'/c 2 ) 

7 * n <v o 2 r h 

u = T + ( 7 nC v , 

■S “V- 

where I is the idem factor dyadic and $ is a unit vector. In our case v is the velocity of the 


comoving system as observed from the laboratory. In the laboratory system an incident plane 
electromagnetic wave is given by: 


E,’ 


( 10 ) 


where x' = x' + a and a ' is an arbitrary constant. Thus in terms of x', ( 10) becomes 


„ , i k • a ' iki'-x'-iWiV 

E.' e ~ e - 1 

<vl 


( 11 ) 


i k'* a ’ 

and e ~ is a constant phase factor which must be carried along. In the comoving system we 
have a plane wave 


„ i k/ • a' iki'X - itojt 
E: e ~ e 


( 12 ) 


where the electric field E ; is determined by the relativistic transformations for the fields. (This is 
not given here because only the phase is of importance. Details are to be found elsewhere 5 ). 
According to the so called principle of the conservation of the phase, the exponents in ( 1 1). ( 1 2) 
must be equal. This prescribes 

k' = U • ( k + to v/c 2 ) k = U * (k' to' v/c 2 ) 


to' = 7 (to + v • k) , to = 7 (to' v • k') . (13) 

We are interested in the effect produced when an incident wave is reflected (e.g., from a plane 
mirror) to a new direction. The reflector is moving according to 

x’ - vf , (14) 

and the local origin x' = x = 0 is chosen such that it is on the plane mirror at f = t = 0. The 
reflected wave is given by 

i k' • u’ i k. , - \ iu)..l 

f () e ~ - e ° . (15) 

where on board the satellite the frequency is unchanged: 
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<*> 0 * Wj , (16) 

and k Q , k| are related by Snell's law; E Q is determined by the pertinent boundary conditions. 
Using the relativistic transformations once more, we obtain in the laboratory system of reference 


E ' 

£o 


. * t i . « # t , #.» 

> kj • » rk 0 *i-iw 0 t 


(17) 


where from (13): 

“o' _ ^ 1 + r£o /c _ 1 ~ l ’ fcj'/ c 

cjj' I kjl 1 + v • kj/c 1 - V * k'/c 

I * <V <N. * *V.V 

Finally, the wave is transformed back to the original system x' , yielding 
H . e ‘< -ho) •« e *JSo -x -iw u ' 


(18) 


(19) 


Eq. (19) is the form taken by an arbitrary plane wave (10) after reflection from a plane surface 
moving at velocity v . 

'■'W 

We now make a few observations. Eq. (18) explains the fact that there is no Doppler fre- 
quency shift in the Fizeau experiment as presented in Fig. 1 . All parts of this system are at rest, 

A 

except the moving fluid. But in the fluid only forward type propagation takes place, i.e.. k 0 = 
a # 

kj . hence according to ( 18) c o Q ' = Wj'. Second, comparing (10) and (19) at x' = 0, t' = Owe 

k ,f — k * s f 

see that the exponential e p ~ describes the round trip phase change due to the initial 

location and the motion of the reflector. Third, since k 0 ' is Doppler shifted with respect to kj' 
this phase factor contains a first order velocity effect. In a velocity independent system. co 0 ’ = cjj', 
hence the time factors in ( 1 0), ( 1 9) are equal and the phase difference between the two waves is a 
constant, independent of time. But when v # 0, then cj 0 ' ¥= Wj' and the ohase difference de- 
pends on time. This is so because the change of phase takes into account the motion of the reflector 
(14) and vanishes only at t' = 0 when the reflector is at x ' = 0. 

In the configuration being considered here, we have waves travelling in both the up- and down- 
wind directions around the triangle. Thus, one must consider the Doppler effect for waves 
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propagating in opposite directions, as depicted in Fig. 2. Assuming S } to be in motion, the Doppler 
effect for rays 1 , 1' is given by (18). For 2', 2" the same formula applies, but the directions and 
frequencies will be denoted by bars. 


A 

“o 1 + x * W c 


“i' 1 + v • Jc,/c 


( 20 ) 


Now align the rays such that in the comoving system they are oppositely directed 


A 

k. 


A 

h 


A 

k. 


A 

k. 


( 21 ) 


Note that due to the relativistic formulas for the aberration phenomenon (13), there will he some 
difference in directions, either in the laboratory or in the comoving system of reference. The pre- 
sent choice simplifies the discussion. Using 


Wj’ = Wj' (22) 

and substituting ( 18), (20). (21). we obtain 
“o' 1 < v • kj/c) 2 

— = . (2.)) 

“o' 1 !io c)2 


Thus, in this case, the reflected rays 1 ” and 2'. bigs. 2,3, do not have the same frequency , the dif- 
ference being of the order ( \ c)‘ • to This is to say that reflection from a moving mirror is non 
reciprocal, in the sense that a net frequency difference occurs when the roles of incident and 
reflected waves is reversed Although this effect is of second order in \ c it can be significant be- 
cause it introduces a time dependence into the problem. Moreovci . it ( 2 I ) is not satisfied, a first 
order effect w ill appear in (23). This can happen as a result of beam spread and orbital motion, 
discussed below 

I- fleets due to orbital motion and beam spread: For a finite re Heel or the above results ma\ he 
still used pronded edge diffraction effects. uv negligible. However, it must be noted that additional 
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geometrical effects are present, which must be adequately taken into account. As the reflector 
moves, it intercepts different parts of the incident beam. Also the outgoing beam is laterally dis- 
placed as depicted in Fig. 4. However, as long as the reflector and the receiver are well within the 
beams, the above plane wave formulas are applicable. Pulses emitted simultaneously from M (in 
opposite directions around the triangle) will reach Sj about 0.01 sec apart because of the different 
paths traversed. During this time the reflector Sj . say, moves a distance of the order of 100 m. This 
distance is small compared to the expected beam cross sections and therefore the plane wave for- 
mulas can be used, provided the directions of the rays are properly taken into consideration. 

Due to the orbital motion of the satellites, the velocity does not remain constant, and the kine- 
matic effects might affect the results of the measurement. For a satellite moving at a height of say 
1000 km the absolute value of the velocity remains practically constant during a period on the order 
of a second. However, the direction of the velocity is changed on the order of 10“ rad/sec. This 
will affect the Doppler shift ( 18) and constitutes a perturbation which must be taken into account 
in the signal processing (discussed in the next section). 

The analysis of the Doppler effect given above is based on the assumption of plane wavefronts. 
Due to the spreading of laser beams, wavefronts become curved. In Fig. 5 it is assumed for sake of 
an illustration, that the wavefronts are spherical, having point a as the center of curvature. Ficti- 
tious rays and the fictitious extension of the reflector are shown in dashed lines. The moving re- 
flector first engages rays 1.1'. later it intercepts 2. reflecting it as 2'. By inspection of the fictitious 
rays 1.1" and 2. 2*' it becomes clear that the only effect on our analysis is again a change of diroo- 

A 

tion. this time for the unit vectors k in ( 18). This effect is of the same magnitude as the above 
orbital motion effect. 

Hffect of relative motion on measurements: In general both satellite and atmosphere (wind) 
are in motion relative to an observer on the earth. Thus the question of what velocity will register 
in our measurements is of importance. To be specific, we first consider the Fi/cau experiment ol 
Fig. 1 , assuming v = 0. The phase accumulated by the raysdependson the electrical length (i.e.. the 
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equivalent free space length) of each path. The phaae difference provides the zero reference for our 
experiment. Now let us set the upper vessel (Fig. 1) in motion, at a relative motion v R with respect 
to the laboratory. Because of the ensuing Doppler effect, the excitation frequency for the upper 
vessel is now different. This is tantamount to saying that the electrical length is modified. This will 
shift the zero reference, but will not otherwise affect the results of the Fizeau experiment. This is 
due to the fact that the Fizeau effect is already of first order in v/e, hence changes of frequency due 
to motional effects are negligible in (7), (8). A detailed .'.iscussion. pertinent to the configuration of 
Fig. 2 is given in the next section. An investigation of relative motion in the Fizeau effect has been 
conducted by Zeeman (e.g.. see Jones 4 for reference to original papers, see also Zernike 9 ). 
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SIGNAL PROCESSING 


In order to extract the phase change due to the motion of the atmosphere, it will be necessary 
to compare the phase of the two laser beams which have propagated around the path (Fig. 2) in 
opposite directions. First, consider the case where the atmosphere is absent. To further simplify 
the analysis, it is temporarily assumed that orbital motion and beam spread effects are absent. These 
restrictions will be waived later on. To begin, consider the phase accumulated by a plane wave 
transmitted from the master station, M, and travelling around the path in the direction M, Sj , S 2 , 

M. (From now on primes wilt be suppressed since only observations in the laboratory frame will 
be considered.) Assume that the wave transmitted from M is: 


£M1 e 


(24) 


where k MJ is directed from the local origin x^ = 0 to the local origin x j = 0 of Sj . Consequently 

i k ui * lug 

the phase factor e " is introduced in translating the wave to the local origin of Sj . Al- 
though the choice of x . = 0 is arbitrary, if we choose it on the actual reflector at t = 0 (and not 
* 

its fictitious extension as in Fig. 5) then k^j and a^j are parallel and k^j * a M1 - k M) a M] . The 
wave reflected from Sj is therefore given by 

E 12 e ikM,aM1 , (25) 


where k jt is directed towards the local origin X 2 = 0 of S 2 . The wave arriving at S 2 will have the 

i k 3 

additional phase factor e 12 12 , where is the distance from Xj = 0tox2 =0. Thus, the 
wave reflected by S 2 in the direction of M is given by 


p ( k M! a Ml + *12*12) 

£2M e 


e > *2M ‘ *2 " »W2M t 


(26) 


where k 2 ^ is directed towards the local origin of M. Hence at M we receive at x^ =0 a signal 


,.(a) >[ k Mt a M , ♦ k, 2 a 12 + *2M a 2M “ W 2M •) _ ,,(») i 

Lli V — Cw e 

.•*^**1 


(27) 


15 



where is its amplitude. Similarly, for the round trip 2, 2', 2" we start with a wave emitted by 

M, 


where by now the notation is obvious. The wave reflected by Sj towards Sj will be 
e ik M2*M2 


(28) 


(29) 


There v iil be a similar reflection produced by Sj , and so finally the wave arriving at Xj^ = 0 is 

Jb) i[kM2*M2 + k 21*21 + k lM*lM“ w lM t ) _ r.( b ) J ^ b , 

£M e = £M e ( 


which will be compared to (27) in order to extract the phase difference due to the Fizeau effect. 

Notice that k ij # k ji (e * ’ k MI * k lM>’ since the Doppler frequency shifts are different in each 

(a) (b) 

case. Also note that the amplitudes £j^| and are unimportant for our problem and so no 
explicit expressions are included for them. Any boundary condition at the reflectors Sj and S-, 

would be identical for the two oppositely traveling waves and is already absorbed into the ampli- 

, < r< a) i r (b > 
tudes and Ey . 

In view oft 18) all and Wjj in (27). (30) are proportional to the source frequency cc. Hence 
^ a . ( 27). and (30) can be recast in the form 
0 a = co (a + 0t) , 

0^ = co (a + 0t) (31) 

Now suppose that the two beams a and b are coherently detected and the phase difference AO 
0 a 0^ measured, yielding: 

AO = co | (cv a) + (/? /?)t|. (32) 

I’aking A “ lO/i. and v s - 10 knt/sec for the satellites, to be representative values, according to 
( 23 ) we tind the difference frequency ujjhi = cj( 0 /J) to be on the order ol cj (v s "c)- ~ 

1 00 k 11/ Next we consider the values of co (a a). ( 32). For an observer on the ground we have. 
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to the first order in the velocity (Fig. 2) : 

k Ml " k < ! + 'hil^ " + °Ml> * 

k M2 * k < 1 + *M ‘W 0 ) * k <* + °M3> > 

k 12 = k Ml “ ll * lLM|/ c + k 12 *Il/ c > " k Ml< ! + ®l)* 

k 2M * k 12< 1 ~l2 * k 12 /c + l2 “ k, 2 (l+o 2 ), 

k 21 " k M2 < ! ~l2 * £m2^ c + *2 * k 2l/ c > “ k M2 <* + °2>* 

k lM = k 21 0 -*1 * k 21 /c + *1 ,k lM/ c > * k 21 (l+o,), (33) 

and for ay = a^j (e.g. , a M j = a JM ) this yields 

<*>(a - a) * k U M , (o M , - <*m 2 - o 2 - o,) + a M2 (o M , + Oj + o 2 - o M2 ) 

+ a 12^°Ml + °l °M2 °2^ - k = w / c - (34) 

Taking ay on the order of 10 4 km, it turns out that w (a - a) is on the order of 2ir • 10 7 rad. 

Since u> (a - S) is much larger than A0 which we are trying to measure, a method must 
be found to calibrate the system for w(a-a) prior to measurement, and to account for 
its variation during the measurement. This will be discussed subsequently. 

There exist a variety of effects changing w (a - a) and u(0-0) from the values obtained 
above. Note that signals simultaneously emanating from M do not reach Sj (or S 2 ) at the same time 
since the paths that they follow are of different length (e.g., MSj Sj compared to MSj). The time 
difference is on the order of 10" 3 to 10 -2 seconds. As a result the o, in (33) associated with k, 2 
and k, M arc not the same, nor are the o 2 associated with k 2 , and k 2M identical. Taking this effect 
into account only slightly changes the magnitude of w (a ~ a) in (34). The most important effects 
are the changes that occur in satellite velocity with time as the satellites move in their orbits around 
the earth (orbital motion), and the fact that the laser beams are not truly plane waves (beam spread). 
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These two effects can introduce a time dependence in Ojh , . Oj . Oj • To study the effect of 

orbital motion, let us represent (34) by 2* * 10 7 cos 8 , where 8 is the angle understood in the 
scalar products in (33). Due to orbital motion 8 is time dependent and therefore a representative 
for (34) can be written in the form 

2ir • 10 7 cos(0 n + r) , (35) 

where r is the distance of the satellite from the earth's center and v t is its velocity, and 0 O refers to 
the time r * 0 which is arbitrarily chosen. Thus during one second v. changes by an angle on the 
order of v s /r » 10“^ rad. For worst case analysis take | 0 Q | = itj 2, hence near r = 0 (35) behaves 
as 2ir • 10 4 t, corresponding to a frequency on the order of 10 kHz. The time dependence intro- 
duced by beam spread is of the same nature and magnitude, with r - 10 4 km in (35) standing for a 

A A 

typical distance between satellites. The effect on w (0 - 0) is even larger, since 0 0 corresponds 

to a small difference between large numbers. Here the time delay for simultaneously emitted signals 
arriving at a reflector introduces a change of direction which perturbs(21). Hence (23) contains 
first order effects which must be carefully evaluated. For v s * 10 km/scc. a distance 1 0 4 km be- 
tween satellites and delay time on the order of 10"^ seconds, the angle subtended by the oppositely 
propagating beams is on the order of 10"^ rad. This number is on the order of v/c hence the fre- 

- A 

quency difference is on the order of (v/cp as in (23); i.e., u)(0 0) is still considered to be in the 

A 

100 kHz band. Depending on the parameters chosen in Fig. 2. o>X0 0) on the order of 1 MHz can 
also be considered realistic. In any case, for a short period of observation, on the order of 1 10 

seconds, it can be assumed that the time dependent effects combine to yield in (32) a fixed value 

A ^ 

oj (a am * a sl'phtly modified frequency u (0 0). Of course, we cannot hope to compute 

these values from the formulas given above with sufficient accuracy allowing for measurement of 

A 

the Fi/cau effect, however u> ( 0 0) is amenable to very high precision measurement . ami lor 

A A 

u) (a °d T r() 0| dy its deviation front an integral number of 2ir is relevant < re., uj (o «) rr Q ,m>l Fdo 
2n) Commercially available counters are capable of time measurement with an error on the order 
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of 10“ 10 , which wit) also be the error in frequency for a one second measurement (atomic clocks 
are several orders of magnitude more accurate). This suggests a method for calibrating the system. 
Thus if we measure the frequency w (0 - 0) and make a phase measurement at some time r = 0 
(prior to the time when the line of sight penetrates the atmosphere), then we can safely assume that 
u) (o - a) rm Q + u> (0 ~ $) r is known for a period on the order of 1 - 10 seconds. 

This brings us to the point where we wish to consider the effects of the atmosphere on our 
system. There are three effects taking place simultaneously, and since the Fizeau phase (7), (8) is 
the one we want to measure, the other two must be computed, using independent data. The first of 
these two effects is the relative motion in the system, which produces different excitation frequen- 
cies for the oppositely moving rays. In (27) we have to subtract the free space phase k^j L and 
add the effect of the refractive medium as k MI Ln. Similarly k JM (n - 1) L is added to (30). 
The difference is given by 

AiJ/r = k ( 0 M1 - o 1M ) LN • I0" 6 , (36) 

which takes into account the relative velocities, as explained in the previous section. In (36) the 
same effective value LN is assumed for the two oppositely going waves. Comparing (36) and (7). it 
is seen that the Fizeau phase is three or four orders of magnitude smaller, because of the ratio of 
the wind velocity to satellite velocity. However, satellite velocity can be accurately measured by 
monitoring the motion. It can also be inferred from the satellite's height, assuming a circular orbit. 
For example, if the positioning error is 10 m, we use v t “ = rg, where g * 10 m/sar is the gravita- 
tional constant relevant to the distance r from the earth's center, obtaining Av s /v s * 10“ 6 . To 

A A 

compute (36) we need the directions kj^j . k j^f . With a positioning error 10 m and distance be- 
tween satellites on the order of 10,000 km, the error is on the order of 10~ 6 . Satellite positioning is 
continuously improved, and with the advent of the Global Positioning Svstem project one or two 
orders of magnitude improvement can be expected. Hence in computing (36) we still have enough 
precision left for determining the Fizeau phase (7). (8). The second effect which must he computed 
independently is the change in path length and direction of propagation of the line of sight uue to 
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ray bending in the atmosphere. In aider to assess *his effect we assume that on a path length of 
100 km in the lower atmosphere the ray bending will be significant. Using the (4/3)a effective 
earth radius model 6 , this means that the radius of curvature of the ray will be 4a. This will bend 
the ray through an angle on the order of 10"’ rad at the extreme level of penetration into the atmo- 
sphere. when the line of sight is close to occultation. The effect on the path length is on the order 
of I cos 10' 3 * I0" 3 hence with w (a - a) * 2* • 10^ we have to compute a phase on the 
order of 20* which is much smaller compared to A^g . On the other hand, the change of direction 

i A 

on the order of 10“ J rad implies a change in the Doppler effect factors k • v/c. The effect on 

A A A 

u) ( 0 0) is negligible, because 0, 0 are affected in the same way, leaving the difference 0 0 prac- 

tically unaltered. However, the factors o in (33) are changed, implying a phase on the order A^g = 
2* * 10 4 which is larger than A^g. It is therefore expected that the minimum altitude for measure- 
ment of wind speed w ill be limited by the accuracy with which Avi'g can be determined. 

Prior to the line of sight penetrating the atmosphere A y is measured and its value during the 
measurement can be projected by knowing the frequency, as explained above. While the remote 
sensing measurement of the wind velocity is taking place, a different phase A^ is recorded. The 
difference, including the correction factors A^g. A^g discussed above finally yields the Fi/.eau 
phase A0. according to 

A# = A^^ (Ai^ + A^g + Avi'g). (37) 

From this the wind velocity is obtained as an average value over the path, for various altitudes. 

We have considered the various effects separately; the combined error in measuring the wind 
velocity is expected to he larger. On the other hand, worst case parameters have been chosen most 
of the tunc, and some effects will tend to mutually cancel A better understanding of the interac- 
tion of the various effects calls for a computer modelling of the entire problem. 
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COHERENCE CONSIDERATIONS 

Since the proposed procedure for measuring wind velocity depends on making phase measure- 
ments, the coherence of the User beams as they traverse the system must be examined. 

The problem will be discussed below in a preliminary way. Clearly, better understanding of this 
aspect of the problem is needed. 

First, the coherence time of the source must be examined. In the configuration of Fig. 2 ray I 
reaches Sj before 2'. Assuming the distance difference to be 3000 km, this amounts to a dcUy time 
of 0.01 sec. During this time S| was moving at a velocity of 10 km/*ec, say, covering a distance of 
100 m. It follows that the time difference between rays 1", 2" returning to M, Fig. 2, is of the 
order of 1 nsec. Such a coherence ' 'me is of course amply supplied by a User. 

Next, the randomness of the a>mosphere must be taken into account. Over the long paths 
pUnned here, it is expected that atmospheric turbulence and irregularities will degrade the phase 
information. What we shall argue here is. that this does not invalidate the fundamental ideas given 
above, and that there are ways of controlling the amount of incoherence encountered at the detec- 
tor. 

In order to understand the physics of the problem, consider a point source in a random me- 
dium like the atmosphere. Due to the random irregularities in the atmosphere, the wave will develop 
phase fluctuations, i.e., the wavefront away from the point source will not be spherical any more. 
The presence of a distorted wavefront also means that rays, perpendicular to wavefronts, will now 
travel in directions which are not strictly radial. At Urger distances these rays will interfere, giving 
rise to the scintillation phenomenon. This problem has been studied extensively, both theoretically 
and experimentally. See Ishimaru^ , and Tatarski 1 ®, who also cite many earlier references. Since 
the present problem is closely related to scintillation from point sources, ideas relevant to this sub- 
ject will be used. 

A discussion well suited for our problem is given by Lawrence 1 1 . He argues that the irregulari- 
ties most effective in producing scintillation are of dimension of the first Fresnel /one. Hence, in 
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order to use point source theory, at least the first Freanei zone for X * I Op must be illuminated. 

For the present parameters where satellites arc about 1000 km high above the ground and about 
7500 km apart, the Fresnel zone is of the otder of 4 m. midway between the satellites. Typical 
laser beam spread angles are of the order of X/D. D being the aperture diameter of the transmitter. 
Taking D to be 0. 1 m. we obtain an angle of I0“* rad. This means that the radius of the beam's 
cross f ^tion will be hundreds of meters, containing many Fresnel zones. We arc therefore justified 
in treating the radiation as originating from a point source. The diameter of the most effective irre- 
gularity, which is also the radius of the first Fresnel zone, is given by 

d = -yAg . g = Zj /.,/ U, +z 2 ) (M 

where /j . are the distance to the source, the distance to the receiver, respectively. As a first 
approximation for our case, the atmosphere -an be thought of as being lumped midway between 
the satellites, so that z ( = /, and as mentioned above, d * 4 m. It is argued that in the presence of 
larger irregularities, the smaller ones predominate, much like the case of a ground glass plate put in 
front of a lens. On the other hand, the spectrum of atmospheric turbulence increases with irregu- 
larity size. It is therefore assumed that irregularities whose size is given by ( 38) are most effective 
It is also known that the validity of the present model is limited * *. C orresponding to the irregular- 
ities a random pattern of intensity fluctuations will be measured in the vicinity of the receiver. I he 
pattern radius corresponding to (38) is given by* * 

p ~ <1 + z, z,ld 2. t.V)i 

If we take Zj = /, again, (39) yields 

p ' d * 4 m (40) 

The pattern radius is closely related io the distance between uncorrelated parte in the held 
measured near the receiver This concept is important tor our discussion of apeture averaging, 
given below. I awrcnce* * displuvs the normalized covariance of log-intensity fluctuations as a 
tuiKtionotf = AtffXL . where l is the tine of sight path length and A is the dislam between 
two small-aperture detectors for out present assumption o! d \ XI wc have 



f - (6/d) ir"* 4 * 6/(2d). 
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It to realized that the theory corresponding to Fig. 25.6* 1 to more intricate, however, we use the 
results here to get a rough idea of the parameters involved. Thus by inspection (Fig. 25.6* 1 ), it is 
found that f ■ I roughly corresponds to zero correlation, i.e., i) ■ 0, and 

6 » 2d * 8m (42) 

between detectors will ensure uncorrelated statistical measurements 

The problem of decreasing scintillation is very similar to overcoming fading in radio wave 
propagation, and one obvious method is spatial diversion reception. The analog for the present 
problem is the increasing of the receiver’s aperture. Loosely speaking, if the aperture collects more 
rays, at different phases, the random phase factors will be eliminated, and the coherent component 
of the radiation will be enhanced. This is usually referred to as aperture averaging and is discussed 
by ishimaru^ and Tatarski 1 ®, for example. For the aperture averaging to be effective, enough un- 
correlated “portions” of radiation must be added through the aperture. In systems where good 
imaging quality is also required, this implies large aperture, costly and heavy telescopes. The imag- 
ing problem docs not enter into our considerations, therefore for the present system a large array of 
small aperture receivers will suffice. This poses the unrealistic requirement of having an aperture 
many times larger than 8 - 8 m. Fortunately, the present configuration corresponds to a detector 
moving at about 10 km/sec. This means that during one second as many as 1 250 uncorrclatcd 
samples can be gathered. Of course, we are limited by the fact that the satellites change position 
continuously, sweeping through different parts of the atmosphere. However, during the time of the 
order of one second, the distance traversed is of the order of a few kilometers. This still allows for 
good resolution of the order of 10 km. The error in phase due to scintillation may be considered 
as noise present in the process of measuring &$. The present method of “synthetic” aperture aver- 
aging will impinve the signal to noise i.itio It y taking the average ol a lew huiulied samples. £<!> will 


be enhanced and the noise* diminished 



\ 


The above is a somewhat pessimistic evaluation of the system, since there arc additional factors 
working in our favour. For example, we have to take into account the fact that the transmitter is 
moving too, continuously changing the line of sight, at different positions and angles. This should 
have an effect at least as significant as the moving receiver. Consequently 6 = 4 m and twice as 
much independent data can be accumulated in a given time period. It is well known that scintilla- 
tion is rapidly reaching saturation, and does not grow with the length of the line of sight 
path 3 .1 0 .1 1 >1 2 . Lawrence 1 1 puts the saturation distance for visible light at a path near the ground 
in the vicinity of one kilometer. This phenomenon, combined with the low scattering in the atmo- 
sphere in the lOg IR band might result in very low noise levels. The subject will have to be discus- 
sed in more detail. 
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SUMMARY 


A method has been proposed for measuring atmospheric wind velocity from space platforms. 
The present method is based on the Fizeau effect, and consists of transmitting two laser beams 
through the atmosphere, one upwind and the other downwind, and measuring the phase difference. 

Typical numbers for the atmosphere indicate that the effect should be measurable. The ques- 
tion of carrying out the measurement in the presence of Doppler effects has been considered in 
detail, and it has been shown that although these effects make the measurement more difficult, the 
wind velocity can be measured, in spite of the fact that it is several orders of magnitude lower than 
the satellite velocity. 

Inasmuch as a coherent measurement is proposed, the mechanisms introducing incoherence 
have been discussed. The main effect is expected from scintillation, which can be decreased by 
averaging the measurement over a short time period. 

The present study is only preliminary, and many questions must still be answered to determine 
the practicality of this method for measuring wind velocity. 
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Figure 1 . Schematic of the Fizeau experiment. Two coherent light beams traverse a moving liquid 
in opposite directions. The phase difference measured at the output is proportional to the 

medium’s velocity. 




Figure 3. Estimation of atmospheric path length. 
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4. Effect of reflector motion on lateral displacement of beams. 
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APPENDIX A: A LIST OF OPEN PROBLEMS 


1 . Detectability of Signals : Existing and projected instrumentation must be examined for imple- 
mentation of the present method. This applies especially to lasers and detectors. 

2. Atmospheric Attenuation : Based on available data for standard atmospheres and special deviat- 
ing cases, estimate the bounds on the applicability of the present method. 

3. Cloud Cover: Using climatological data, estimate the fraction of time for which the method is 
operable. Compare to limitations on other methods. 

4. Oce an Reflection Configuration : This configuration has been abandoned in the present study, 
because it was felt it might be too noisy. The method is attractive because it requires two satellites 
only. Check if aperture averaging technique, as described above, facilitates the use of this configu- 
ration. 

5. Coherence Problems: Provide a quantitative analysis of the incoherence introduced bv the 
atmosphere and the reduction of scintillation by synthetic aperture averaging. 

6. Resolution: C onsider inversion techniques relevant to the present system for improving resolu- 
tion. 

7 . Compute' Model: The combined effect of many factors described above is too complicated to 
be investigated unalyt ieallv. A computer model should be constructed in order to test various ideas 
given above. 

8. An Acoustical Analog for Doppler Effects: The design of a laboratory or ground based experi- 
ment w hich car. simulate the high velocities of the satellites is probably as complicated as using the 
system itself. In acoustics it is relatively easy to achieve strong Doppler effects (although we cannot 
simulate the second order relativistic effects discussed above). Design an appropriate experiment, m 
air or water, that will test the feasibility of coherent measurements in the presence of strong Dop- 
pler cl feels 
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